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SUMMARY 

Computational movies were used to analyze the formation and collapse of 
vapor cavitation bubbles In a submerged journal bearing. The effect of 
vibration amplitude on vapor cavitation was studied for a journal undergoing 
circular whirl. The boundary conditions were Implemented using Elrod's algo- 
rithm, which conserves mass flow through the cavitation bubble as well as 
through the oil-film region of the bearing. The vibration amplitudes for the 
different cases studied resulted In maximum eccentricity ratios ranging from 
0.4 to 0.9. The minimum eccentricity ratio reached In each case was 0.1. For 
the least vibration amplitude studied In which the eccentricity ratio varied 
between 0.1 and 0.4, no vapor cavitation occurred. The largest vibration 
amplitude (l.e., eccentricity ratios of 0.1 to 0.9) resulted In vapor cavita- 
tion present 76 percent of one complete orbit. 


INTRODUCTION 

Vapor cavitation In submerged journal bearings Is a type of cavitation 
that is only present when the bearing Is In a dynamic state. Normally vapor 
cavitation appears during some part of a vibration when the dynamic loading on 
a bearing results In a rapid Increase of the minimum oil-film thickness. Sur- 
face separation due to normal motion causes film pressure to decrease until 
the tensile strength of the oil or the bonding strength of the oil to the 
bearing surfaces Is exceeded. The cavitation bubble produced contains mostly 
oil vapor of a very low pressure (close to vacuum). 

Quite frequently, vapor cavitation produces erosion on the bearing sur- 
faces when the vapor bubbles are transported to a high-pressure area of the 
lubricating film where they collapse. The presence of vapor cavitation also 
changes the dynamic behavior of the bearing. In particular, the damping In 
the bearing decreases considerably. 

The phenomenon of vapor cavitation In journal bearings was experimentally 
studied by Jacobson and Hamrock (refs. 1 and 2). They found that the size of 
the vapor cavitation zone was dependent on factors such as the vibration 
amplitude, the minimum oil-film thickness reached during a vibration, and the 
surface energy of the bearing surfaces. 



This report Is a theoretical-numerical Investigation of the appearance, 
growth, collapse, and disappearance of the vapor cavitation zone for a sub- 
merged Journal bearing. A computer program was developed (ref. 3) which uses 
Elrod's algorithm for the continuity of mass flow Into, through, and out of the 
vapor cavitation zone. This program makes It possible to predict the dynamic 
behavior of the vapor cavitation zone. The Cray computer at NASA Lewis 
Research Center was used to do the computations. 


SYMBOLS 

0 shaft diameter, m 

e eccentricity 

g switch function (cavitation Index) 

h film thickness, m 

L bearing length, m 

L/D length to diameter ratio 

m lineal mass flux, kg/m-s 

P fluid pressure, N/m 2 

P a ambient fluid pressure, N/m 2 

P c cavitation fluid pressure, N/m 2 

R radius of shaft, m 

&R radial clearance, m 

t time, s 

U sum of surface velocities In x-dlrectlon, m/s 
V sum of surface velocity vectors, m/s 

W load capacity, N 

w squeeze-velocity, m/s 

x coordinate along circumference, m 

y coordinate normal to x,z plane, m 

z axial coordinate, m 

0 liquid bulk modulus, N/m 2 

y angular position of minimum film, rad 

c eccentricity ratio, e/AR 
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6 fractional film content In cavitation zone; density ratio (p/p c ) In 
full-film zone 

v dynamic viscosity, N-s/m 2 

p fluid density within full-film zone, kg/m 3 

p c fluid density within cavltated zone, kg/m 3 

<p angular coordinate relative to minimum film line, rad 

ud orbital angular velocity of journal center about fixed point relative to 

housing center, rad/s 

u s angular velocity of Journal about Its own center, rad/s 

Subscripts: 

max maximum 

min minimum 


BACKGROUND THEORY 


Elrod and Adams (ref. 4), Implemented a moving boundary scheme that avoids 
Interface complications. The following Is a review of some of the underlying 
factors that led to the algorithm which was later modified by Elrod (ref. 5) 
and used here. The conservation of mass can be written as 

Hf 61 ♦*•(«).<> (i) 

where m represents the lineal mass flux and Is given by 


j» hV _ fih! ft, 

p 2 12y 


( 2 ) 


Substituting equation (2) Into equation (1) yields the Reynolds lubrication 
equation. This equation Is made applicable to the cavitation region as well as 
to the full-film region by Incorporating a switching function that automat- 
ically satisfies the boundary conditions at a moving Interface. Furthermore, 
the fractional film content 0 Is made the dependent variable. This requires 
giving e a dual Interpretation. That Is, In the full-film region, e repre- 
sents the mass content of the film that exceeds the content at the cavitation 
pressure P c . In other words. 


( 3 ) 

where p c Is the density of the liquid at P c . Furthermore, P and p are 
related through the equation for the liquid bulk modulus according to 
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( 4 ) 



or 


P * P c + 8 In e e > 1.0 (5) 

In the cavitation region e < 1.0, e determines the mass content p c eh, which 
can manifest Itself In the form of a smeared mass or striated flow extending to 
both surfaces In the film gap. The fluid transported through the cavitation 
region In the form of an adhered film can be neglected for heavily loaded 
conditions. 


A universal differential equation Is made possible by linking the solu- 
tions of the full-film region with the solutions In the cavltated region by a 
single dependent variable e. However, a cavitation Index, or switch func- 
tion g, was Included so that the resulting partial differential equation would 
be consistent with the uniform pressure assumption within the cavltated region. 
The switch function Is defined from a knowledge of e. Thus, 


g = 0 e < 1 .0 

g = 1 e > 1 .0 


(6) 


The switch function g Is made a factor of the pressure gradient term In 
equation (2) so that the flow Is strictly Couette In the cavltated region. 
Expressing the lineal mass flux m In terms of e and g and substituting It 
Into equation 1 yield the universal differential equation obtained by Elrod and 
Adams (ref. 4): 

♦ \ . [^(eh) ] = 3 Bg(e)fo] (7) 

In the full-film region the solutions of e together with equation (5) 
determine the pressures. In the cavitation region, g = 0 and equation (7) 
becomes 

HF * ! • ■ o < 8 > 


Equation (8) governs the transport of the fluid through the cavitation region. 

The numerical details are given In Brewe's work (ref. 3) as well as In 
the original publications by Elrod and Adams (ref. 4) and Elrod (ref. 5). The 
method combines a control volume approach to deriving the finite-difference 
equations with the use of a switching function to automatically satisfy the 
boundary conditions at a moving Interface. An alternating direction Implicit 
scheme Is applied. At each half-time step the Implicit Euler method Is used 
to march In time. This method Is known to be unconditionally stable (ref. 6) 
Insofar as the ordinary differential equation Is stable. In the first half- 
time step, spaclal differencing along the circumferential direction results In 
a periodic trldlagonal matrix. The solutions for e and, consequently, p are 
found by performing a Gauss-Jordan reduction on the matrix along with a maximum- 
pivot strategy to reduce the error. Spaclal differencing along the axial direc- 
tion Is performed to complete the time step. The solutions for this sweep are 
more readily obtained by using a nonplvoted Gaussian elimination procedure (trl- 
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diagonal solver). Since the journal and housing are aligned, the pressures 
must be symmetric about the axial center. Consequently, the calculations are 
only made over half the axial length of the housing. 


RESULTS ANO DISCUSSION 
Cavitation and Pressure Distribution 

The effect of vibration amplitude on vapor cavitation was studied for a 
journal undergoing circular whirl. The journal bearing configuration Is shown 
In figure 1. The range of amplitudes and conditions considered In this study 
are listed In table I. 

Figures 2 to 4 Illustrate the pressure history within the fluid-film 
region for a journal completing one orbit of Its cycle. Figure 2 represents 
the case In which 0.1 < c < 0.9. The maximum vibration amplitudes In figures 3 
and 4 resulted In e max of 0.7 and 0.5, respectively. As the maximum 
vibration amplitude Is decreased (l.e., c max = 0.9, 0.7, and 0.5) the size 

and life of the bubble decrease (l.e., t = 49.4 ms, 42.1 ms, and 20.5 ms) 

respectively. For e max = 0.4, cavitation was not present. 

Figure 2(a) represents the position of the journal within the housing and 
the associated pressure distribution at the Instant that the eccentricity Is a 
minimum. The dark oblong region Indicates the shape and position of a vapor 
bubble downstream of the minimum-film line. The position of the bubble within 
the small clearance space of the journal bearing Is also Indicated. At this 
particular Instant the tendency of the pressure flow along the side of the 
bubble to cause the bubble to collapse Is becoming greater since It Is propor- 
tional to h 3 . The journal center moves In a clockwise orbit from figure 2(a) 
to (b) (nearly one-quarter of a revolution). While the eccentricity Is 
Increasing, the minimum film Is decreasing - giving rise to the squeeze action 

In the vicinity of the minimum-film line. A slight rise In the pressure field 

can be observed there. Note however that the vapor bubble has moved further 
downstream from Its position In 2(a). It Is now located at the position of 
maximum film thickness. The Increased side flow will cause total bubble 
collapse. 

From figure 2(b) to (c) cavitation Is absent. In figure 2(c), a new cavi- 
tation bubble appears and the pressure generation In the region of the minimum- 
film line becomes significantly large because of the squeeze and wedge effect 
taking place there. Figure 2(d) represents the configuration of the bearing 
when the maximum pressure occurs. Note that this happens 3.5 ms before the 
maximum eccentricity Is reached because of the motion of the shaft (fig. 5). 

As the journal approaches the end of the squeeze cycle and begins to separate 
from the shaft, the sliding velocity approaches a minimum value and the squeeze 
velocity passes through zero. The effect this has on the pressures more than 
offsets any Increase In pressure that would have been realized because of the 
smaller minimum film thickness. 

Figure 2(e) represents the end of the squeeze cycle Illustrating this 
reverse effect. Now the minimum film thickness begins to Increase, creating a 
suction effect, an accompanying growth In the bubble, and a complete dissipa- 
tion of the pressure hump noted In 2(f). It Is Interesting to note that the 
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eccentricity Is the same (0.83) for figures 2(c) and (f). The drastic dif- 
ference In the pressure plot Is a result of the squeeze velocity having 
opposite signs In the two situations. The minimum film thickness continues to 
Increase as It moves In a clockwise direction. The representation shown In 
figure 2(a) Is repeated as Is each of the representations for each succeeding 
orbit. 


Load Capacity Versus Eccentricity 

Figures 6(a) to (d) represent the load capacity as a function of eccen- 
tricity for different vibration amplitudes. The predominant effect to note In 
these figures Is that the cavitation produces a hysteresis effect In the load 
cycle. As previously noted, the size and life of the bubble decrease with 
decreasing vibration amplitude. The hysteresis effect Is more pronounced for 
higher vibration amplitudes. In figure 6(d) the maximum eccentricity reached 
In the load cycle was 0.4 and cavitation was absent. The absence of cavita- 
tion resulted In a load cycle without hysteresis. 


CONCLUSIONS 

The following conclusions can be stated for those cases In which the 
vibration amplitude Is Increasing: 

1. The vapor bubble life Increases. 

2. Maximum vapor bubble size Increases. 

3. The bubble Inception comes earlier and the bubble collapse comes later. 

4. The variation of the load-carrying capacity Increases. 

5. The load cycle changes from nonhysteretlc to hysteretlc when vapor 
cavitation occurs. 
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TABLE I. - OPERATING CONDITIONS 


Radial clearance, AR, rrnn 0.5 

Radius of shaft, R, mm 42.5 

Length to diameter ratio 1.0 

Eccentricity ratio, e 0.1 to G max 

Maximum eccentricity ratio, e max 0.4 to 0.9 

Angular velocity of journal about 

its own center, a> s , rad/s -19.5 

Orbital angular velocity of journal 
center about fixed point relative 

to housing center, wd, rad/s -92.7 

Liquid bulk modulus, B, N/m 2 1.72x10^ 

Dynamic viscosity, y, N-s/m 2 0.066 

Ambient pressure, P a , N/m 2 1 . 01 33x1 0 5 

Cavitation pressure, P c , N/m 2 0.0 




y 

Figure 1.- Journal bearing configuration. 






Figure 3 .- Pressure distribution and bearing configuration for a full period of shaft whirl, ( 0 . 1 ^ £ < 0 . 7 ). Figures (a) to (f) viewed clockwise are consecutive in 
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Figure 5.- Surface velocity sum (U) and squeeze ve- 
locity (w) during a full period of shaft whirl. 



Load capacity. 
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